Effect of macrofaunal bioturbation on bacterial distribution in marine sandy sediments, with special reference to sulphur-oxidising bacteria

We have studied the impact of the bioturbating macrofauna, in particular the lugworm Arenicola marina and the bivalve Cerastoderma edule, on abundances and distribution patterns of total bacteria and of bacteria of selected functional groups in sandy intertidal sediments. The selected groups comprised the colourless sulphur-oxidising bacteria and the anoxygenic phototrophic bacteria, which are expected to occupy small zones at the oxygen–sulphide interface in stable (non-bioturbated) sediments. The presence of a wooden wreck buried in the sediment at 10 cm depth within a large area of intertidal sand flat colonised by lugworms provided a unique opportunity to confront field observations with laboratory simulations. The site with the wooden wreck, which was used as control site, was devoid of both A. marina and C. edule, while the composition of the rest of the zoobenthic community was rather similar to that of the surrounding area. In the field, the density of total bacteria was approximately one order of magnitude higher in the control site than in the natural (bioturbated) site. This can be explained by the higher contents of silt and clay particles (higher surface-area/volume ratio) and higher total organic-carbon contents found at the control site. It appears that the presence of macrofauna affects sedimentation processes, which indirectly influence bacterial dynamics. Samples from the control site have been incubated in the laboratory with A. marina and C. edule added (bioturbated core), while an unamended core served as a control. The laboratory experiments contrasted with the field observations, because it was found that total bacteria were actually higher in the deeper layers of the bioturbated core. Moreover, the populations were more homogeneous (less stratified) and colourless sulphur bacteria were on average less numerous in the bioturbated core. In general, laboratory incubations resulted in a decrease of total bacteria with a concomitant increase of colourless and phototrophic sulphur-oxidising bacteria and thus in modifications of the bacterial community structure. Hence, our results demonstrate that care must be taken in extrapolating results from laboratory experiments (e.g. mesocosm research) to field situations.     

Low densities of sea urchins influence the structure of algal assemblages in the western Mediterranean

Numerous studies of interactions between urchins and algae in temperate areas have shown an important structuring effect of sea urchin populations. These studies focused almost wholly on the effect of high urchin densities on laminarian forests. In contrast, algal communities below 5–6 m depth in the northwestern Mediterranean are characterised by low sea urchin densities (<5 ind m⁻²) and the absence of laminarian forests. No previous research has addressed sea urchin/algal interactions in this type of community. To determine the effect of the most abundant echinoid species, Paracentrotus lividus, on well-established algal communities in this area, we performed a removal–reintroduction experiment in rocky patches located between 13 and 16 m depth in the northwestern Mediterranean, where sea urchin densities ranged between 0.9 and 3.4 ind m⁻². After 6 months, the cover of non-crustose algae was significantly higher in the plots from which sea urchins had been removed than in control plots (84 vs 67% cover). These removal plots reverted to their original state upon reintroduction of sea urchins. The non-crustose algae consisted of turfing and frondose forms, with the former representing some 70% of the non-crustose algal cover. Change in the cover of turfing algae was responsible for the significant increase in algal development in the sea urchin removal plots. The response of frondose algae to the treatment varied between algal species. It is concluded that grazing by P. lividus exerts a significant effect on habitat structure, even in communities with low sea urchin densities, such as those found in vast areas of the Mediterranean sublittoral.     

Simulation of intrinsic random fields of order <Emphasis Type="Italic">k</Emphasis> with a continuous spectral algorithm

Intrinsic random fields of order k, defined as random fields whose high-order increments (generalized increments of order k) are second-order stationary, are used in spatial statistics to model regionalized variables exhibiting spatial trends, a feature that is common in earth and environmental sciences applications. A continuous spectral algorithm is proposed to simulate such random fields in a d-dimensional Euclidean space, with given generalized covariance structure and with Gaussian generalized increments of order k. The only condition needed to run the algorithm is to know the spectral measure associated with the generalized covariance function (case of a scalar random field) or with the matrix of generalized direct and cross-covariances (case of a vector random field). The algorithm is applied to synthetic examples to simulate intrinsic random fields with power generalized direct and cross-covariances, as well as an intrinsic random field with power and spline generalized direct covariances and Matérn generalized cross-covariance.     

Admissible nested covariance models over spheres cross time

Nested covariance models, defined as linear combinations of basic covariance functions, are very popular in many branches of applied statistics, and in particular in geostatistics. A notorious limit of nested models is that the constants in the linear combination are bound to be nonnegative in order to preserve positive definiteness (admissibility). This paper studies nested models on d-dimensional spheres and spheres cross time. We show the exact interval of admissibility for the constants involved in the linear combinations. In particular, we show that at least one constant can be negative. One of the implications is that one can obtain a nested model attaining negative correlations. We provide characterization theorems for arbitrary linear combinations as well as for nonconvex combinations involving two covariance functions. We illustrate our findings through several examples involving nonconvex combinations of well-known parametric families of covariance functions.     

Assessing the accuracy of sequential gaussian simulation through statistical testing

Sequential Gaussian simulation is one of the most widespread algorithms for simulating regionalized variables in the earth sciences. Simplicity and flexibility of this algorithm are the most important reasons that make it popular, but its implementation is highly dependent on a screen effect approximation that allows users to use a moving neighborhood instead of a unique neighborhood. Because of this, the size of the moving neighborhood the number of conditioning data and the size of variogram range are important in the simulation process and should be chosen carefully. In this work, different synthetic and real case studies are presented to show the effect of the neighborhood size the number of conditioning data and the size of variogram range on the simulation result, with respect to the reproduction of the model first and second-order parameters. Results indicate that, in both conditional and non-conditional simulation cases, using a neighborhood with <50 conditioning data may lead to an inaccurate reproduction of the model statistics, and some cases require considering more than 200 conditioning data. It also can be understood from the result of example 3 that when the variogram range is beg compared to the simulation domain determination of inaccurate simulation program is harder.     

Western boundary upwelling dynamics off Oman

Despite its climatic and ecosystemic significance, the coastal upwelling that takes place off Oman is not well understood. A primitive-equation, regional model forced by climatological wind stress is used to investigate its dynamics and to compare it with the better-known Eastern Boundary Upwellings (EBUs). The solution compares favorably with existing observations, simulating well the seasonal cycles of thermal structure, surface circulation (mean and turbulent), and sea-surface temperature (SST). There is a 1.5-month lag between the maximum of the upwelling-favorable wind-stress-curl forcing and the oceanic response (minima in sea-surface height and SST), which we attribute to onshore-propagating Rossby waves. A southwestward-flowing undercurrent (opposite to the direction of the near-surface flow) is also simulated with a core depth of 1000 m, much deeper than found in EBUs (150–200 m). An EKE budget reveals that, in contrast to EBUs, the upwelling jet is more prone to barotropic than baroclinic instability and the contribution of locally-generated instabilities to EKE is higher by an order of magnitude. Advection and redistribution of EKE by standing mesoscale features also play a significant role in EKE budget.     

Late Quaternary biotic and abiotic controls on long-term sediment flux in a northern Australian tropical river system

The modern distribution of monsoonal rainforest in the Australian tropics is patchy and is mainly associated with river corridors and groundwater springs, which indicates a strong dependence on hydrologic and geomorphic conditions. While their present distribution is well known, very little data exists on past spatial and temporal dynamics of these ecosystems, or their medium- to longer-term controls. Factors such as (i) fire frequency and type, and/or (ii) hydroclimatic conditions (e.g. droughts) have been proposed to control riverine corridor rainforest extent. Recent observations, however, also suggest an additional (iii) geomorphic control induced by alluvial knickpoint migration. Sediment sequences provide valuable archives for the reconstruction of longer-term (a) floodplain sedimentary dynamics, (b) local vegetation history, and (c) catchment-wide fire histories. This study investigates such a sediment sequence at Wangi Creek, and shows that a phase of aggradation, lasting ~4000 years, was recently disrupted by channel incision and floodplain erosion. The aggradational phase is characterized by sand deposition with average vertical floodplain accretion rates of 0.8 cm/yr and includes phases of soil development. The recent incisional phase has changed hydro-geomorphic conditions and caused widespread degradation of vegetation, erosion and lowering of the macro-channel surface. While there is no evidence in our data for an erosional event of similar magnitude since the onset of late Holocene floodplain aggradation, Wangi Creek experienced significant erosion and incision immediately before ~4000 years, providing the first evidence for a tropical cut-and-fill river system. We hence argue that phases of aggradation mainly controlled by biotic processes alternate and depend on feedbacks with incision phases controlled mainly by abiotic processes. The results show that eco-hydro-geomorphic feedbacks may play a crucial role in the medium- to longer-term history of tropical fluvial systems and need to be considered when interpreting fluvial archives with regards to climate, fire or human induced change. © 2019 John Wiley & Sons, Ltd.